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Using high-resolution photoemission spectroscopy we demonstrate that the electronic structure 
of several organic monolayer systems, in particular 1,4,5,8-naphthalene tetracarboxylic dianhydride 
and Copper-phtalocyanine on Ag(lll), is characterized by a peculiar excitation feature right at the 
Fermi level. This feature displays a strong temperature dependence and is immediatly connected to 
the binding energy of the molecular states, determined by the coupling between the molecule and 
the substrate. At low temperatures, the line-width of this feature, appearing on top of the partly 
occupied lowest unoccupied molecular orbital of the free molecule, amounts to only ~ 25 meV, 
representing an unusually small energy scale for electronic excitations in these systems. We discuss 
possible origins, related e.g. to many-body excitations in the organic-metal adsorbate system, in 
particular a generalized Kondo scenario based on the single impurity Anderson model. 

PACS numbers: 73.20.At 
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For more than twenty years there has been a thor- 
ough investigation on 7r-conjugated organic molecules, 
which have shown to be suitable for the application in 
organic electronic devices [I]-^. These molecules often 
form long-range ordered films on single-crystalline metal 
substrates, allowing a systematic and fundamental study 
by various surface sensitive techniques. In particular 
from electron spectroscopy methods, as photoemission 
spectroscopy (PES), inverse photoemission (IPES), and 
x-ray absorption (XAS), deep insight into many impor- 
tant features of the electronic properties of condensed 
films and their interfaces has been achieved [^{5]. The 
latter are of particular importance since they crucially 
determine the properties of possible devices. Additional 
microscopic and spectroscopic information has been ob- 
tained by use of scanning tunneling microscopy (STM), 
leading to complementary information about the relation 
between geometrical and electronic structure [5Hl2j. 

Usually, even the spectra of highly ordered films show 
features with a line-width of several hundreds of meV, 
mostly determined by vibronic excitations within the ad- 
sorbed molecule [71 [T31 [Tl] . Therefore, even experiments 
with modern high-resolution photoemission spectrome- 
ters for VUV photoemission (UPS) display only features 
which are about two orders of magnitude larger than the 
most narrow peaks in other solid state or surface sys- 
tems [T3]. Local spectroscopic measurements by STM, 
on the other hand, show narrow features in the tunnel- 
ing conductivity measurements through a single molecule 
PUI [TBH19) . which have been attributed to a possible 
Kondo like process within the charge transport through 
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the adsorbed molecule. However, neither a direct ev- 
idence of a local magnetic moment, necessary for the 
Kondo effect, nor an immediate experimental observa- 
tion of the spectral function does exist yet. 

Here we report about a high-resolution photoemis- 
sion study on two different organic monolayer systems 
that display a new narrow peak in the excitation spec- 
tra near the Fermi level, possibly related to strong 
electronic correlations in the system. Since the pho- 
toemission spectrum can be interpreted as the single 
particle spectral function A$(E) [T3 [50], it allows in 
general a quantitative determination of the influence 
of many-body effects in the system. After describ- 
ing the experimental setup and the sample preparation, 
we discuss 1,4,5,8-naphthalene tetracarboxylic dianhy- 
dride (NTCDA) on Ag(lll), and compare it to Copper- 
phtalocyanine (CuPc) then, which has an analogous spec- 
tral behaviour, although the geometrical structure of the 
overlayer is different. 

The experiments have been performed on a UHV setup 
based on a high-resolution photoemission analyser (Gam- 
niadata R4000) in combination with a monochromatized, 
microwave driven VUV source (He I: /izy = 21.22 eV, He H: 
40.8 eV). The sample temperature during the measure- 
ments was varied between room temperature (RT) and 
approximately 20 K (see Ref . I^TI for details). Substrate 
preparation (sputtering and annealing j22)). surface char- 
acterization by low-energy electron diffraction (LEED) 
and x-ray photoemission spectroscopy (XPS), and depo- 
sition of the molecules from a Knudsen cell [23] have been 
done at RT and in situ. Due to radiation damage, the 
spectra have shown a decrease of the LUMO (lowest un- 
occupied molecular orbital of the free molecule) intensity 
by about 20% after an exposition time of « 30 min to the 
VUV light. Therefore, we have minimized the exposition 
time to less than 5 min per spectrum, repeatedly reposi- 
tioned the sample to non-irradiated areas, and carefully 
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FIG. 1: (Color online) Series of subsequent high-resolution 
photoemission spectra (Helc, emission angle 35° ) on the re- 
laxed monolayer of NTCDA/Ag(lll) for different tempera- 
tures. The cooling was performed slow (over a period of 3 h) 
leading to an order-disorder phase transition [21]. The top 
three spectra were recorded subsequently at a constant tem- 
perature of 130 K. The solid black line at the bottom indicates 
the contribution of the substrate background. 



checked the reproducibility of the spectra. 

The phase diagram of NTCDA/Ag(lll) shows several 
phases in dependence on temperature and coverage. Here 
we restrict ourselves on the so-called relaxed phase, which 
is characterized by a commensurate superstructure ex- 
plained in detail in [M^EB] representing a net coverage 
of 8 = 0.7 with respect to a dense monolayer [Mj. The 
preparation produces highly ordered monolayers, leading 
to clear LEED patterns and, consequently, to a well de- 
fined backfolding of the Ag 5sp bulk bands in the ARUPS 
data ^27^,28,. 

It was already demonstrated elsewhere [29] that the 
occupied valence regime of the monolayer system is dom- 
inated by the two molecular states HOMO (highest oc- 
cupied molecular orbital) and LUMO or hybridization 
state, in reference to the electronic orbitals of the free 
NTCDA molecule. These states appear above the "back- 
ground" of the metal substrate states, and show a char- 
acteristic angle dependence of the photoemission intensi- 
ties [55] (see supplement), whereas an energy dispersion 
in dependence on fcy can not be observed [2^]. As in the 
related system PTCDA/Ag(lll) HHM ISO] the LUMO, 
unoccupied in case of the isolated molecule, is cut by the 
Fermi edge, demonstrating a partial occupation of this 
state by a transfer of electrons from the substrate to the 
molecule. This is equivalent with a significant hybridiza- 
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FIG. 2: (Color onhne) Analysis of the PES hnewidth 
(FWHM) of the narrow peak at Ef in case of 
NTCDA/ Ag(lll). The spectra were divided by the FDD, 
background corrected and fitted by a Lorentzian afterwards 
(see examples in the inset, normalized to maximum height). 
The resulting line widths are given vs. the sample temper- 
ature T. Assuming an electron-phonon broadening mecha- 
nismus, the linear behavior can be described by a electron- 
phonon coupling parameter of A = 0.68 (black line). The red 
(dark gray) curve was derive from a NCA modelling with 
Tk = 80 K [35]. 



tion of the LUMO with occupied states of the metal sub- 
strate. Note that the LUMO does only appear in the 
photoemission spectra of the molecules in the first layer. 
For molecules farther away from the interface the LUMO 
remains unoccupied [25]. 

The strong interaction between the chemisorbed 
molecules and the substrate leads to a large shift of the 
Shockley-type surface state, which is characteristic for 
the (111) faces of noble metals [22], to energies above 
the Fermi level [HHSS]- In case of NTCDA/ Ag(lll) the 
Shockley state is found 400 meV above Ep [34] . 

Fig. [l] shows a series of high-resolution, normal emis- 
sion PES spectra of 1 ML NTCDA/Ag(lll) for differ- 
ent temperatures. The intensity maxima of the two 
molecular features closest to the Fermi level, HOMO 
and LUMO, appear at approximately 2.4 eV and 0.4 eV, 
respectively. The typical line-width of these molecular 
states is of the order of 500 meV in monolayer systems. 
The background of the substrate is flat near the Fermi 
level and rises steeply at a binding energy of > 4 eV. Al- 
ready in the room temperature (RT) spectrum, the most 
peculiar feature is an additional narrow peak right at 
the Fermi level, appearing on top of the comparatively 
broad LUMO. This new peak can be observed for all 
emission angles and shows the same angle dependence in 
its intensity as the LUMO (see supplement), which next 
to the energy shift described above [M] is another clear 
evidence that it is not related to the Shockley-state ap- 



pearing around the F-points of the reconstructed surface 
Brillouin zones |36j but a molecular feature. Further- 
more, it can not be explained by the molecular states of 
the free NTCDA molecules alone. Therefore, it is rea- 
sonable to suppose that this state is a consequence of the 
close interplay between Bloch states of the substrate and 
molecular states at the interface. 

Very important for the understanding of the new nar- 
row feature is the analysis of its temperature depen- 
dence. Fig. IT] shows a series of normal emission spec- 
tra for temperatures between room temperature (bot- 
tom) and T = 130 K (top three spectra). The sample 
was cooled slowly in this case, i.e. the time between first 
(RT) and last spectrum (130 K) amounts to 3 h. Note 
that the influence of contamination can be ruled out by 
core level PES [53] and X-ray absorption data [53]. The 
spectra display several characteristic changes with de- 
creasing temperature: 1.) HOMO and LUMO shift to 
higher binding energies by about 130 meV, 2.) the new 
peak at the Fermi level becomes narrower and 3.) looses 
rapidly in intensity at temperatures below 170 K until 
the peak has nearly completely vanished at 130 K. Keep- 
ing the temperature at constant 130 K, one can observe 
further changes with the same trend over several hours. 
Going back to RT, the initial spectrum will be restored 
again (not shown). 

The described temperature dependence of the spec- 
tra is immediately related to an order-disorder transi- 
tion [5S]. Below 180 K the long-range order within 
the NTCDA monolayer is distroyed by a thermally ac- 
tivated rearrangement of the molecules. The stable low- 
temperature phase is amorphous, which can be seen by 
the disappearance of both the LEED spots and the back- 
folding of the substrate bands in the photoemission data. 
However, since this order-disorder transition is connected 
to a thermally activated rearrangement of the molecules, 
the transition can be prevented if one cools down rapidly, 
leading to a "frozen" phase with the long-range ordered 
superstructure of the relaxed phase at RT. This frozen 
phase can be stabilized for several hours, which allows 
for a detailed study of the intrinsic line- width of the nar- 
row peak in the ordered phase even down to low temper- 
atures. 

Fig. |2] shows the result of a quantitative line width 
analysis. The individual values for the full width at half 
maximum (FWHM) were obtained by a two step analysis: 
To restore the photoemission signal above Ep, the spec- 
tra have been normalized to the Fermi-Dirac distribution 
(FDD) [m [57PD] first, and then fitted by a Lorentzian, 
which describes the line shape of the peak reasonably 
well (see inset of Fig. [2|. Within the experimental er- 
rors, the narrow peaks appear exactly symmetrically to 
Ep, independent from the temperature. 

The temperature dependence of the line width follows 
roughly a linear behavior, down to minimum values of 
about 25 meV at the lowest temperatures accessible in 
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FIG. 3: (Color online) He I^ spectra (emission angle 50°) of 
Cu-phthalocyanine (CuPc) for two different coverages O = 
0.5 (green (light gray)) and O = 0.8 (blue (dark gray)) at 
T — 80 K. For low coverages, there is no long range order 
formed, as displayed by the LEED patterns {Eo — 12 eV). 
However, there is a resonance peak at Ep. 



our experiment. This value contains the energy reso- 
lution and contributions from extrinsic broadening ef- 
fects as discussed in [22]. Such a linear temperature 
dependence is known from the high-temperature limit 
of the Debye model, where the slope is immediately re- 
lated to the electron-phonon coupling parameter A by 
dV/dT = 2TTXkB [35] ■ From a linear least-squares fit 
of the temperature dependence we obtain here A ~ 0.7, 
which would indicate a very strong electron-phonon cou- 
pling. This is not far from characteristic values of typ- 
ical superconducting metals (A « 1) [H] or C60-systems 
[^2HS] where also a characteristic satellite structure was 
observed. The coupling parameter within the metal sub- 
strate, on the other hand, is only of the order of Aw 0.1 
[15] . Particularly for the similar PTCDA/Ag(lll), elec- 
tron energy loss spectroscopy [46^ indicates substantially 
smaller coupling values than what one can derive from 
the slope in fig. [2] Obviously a strong electron-phonon 
coupling alone can not describe the observed peculiar 
spectral feature in our photoemission data. Note that 
the temperature dependence of the intensity would bear 
additional information but can not be evaluated in the 
present case since the peak is directly at Ep and the nor- 
malisation to the FDD, which is crucial for the analysis 
of the line shape, can lead to artefacts in intensity ^40, . 
To clarify the importance of the long-range lateral or- 
der, we compare the results of NTCDA/Ag(lll) to the 
data of CuPc submonolayers on Ag(lll), which show 
many similarities. At constant temperature this sys- 
tem forms long-range ordered phases for high coverages, 
whereas the individual molecules remain separated and 
disordered for low coverages [IT] HH]. The respective 
spectra at temperatures of T = 80 K are shown in 
Fig. [3] For O = 0.5 (green spectrum on bottom) the 



LEED-picture in the inset shows the disordered g-phase, 
while for 9 = 0.8 (blue spectrum on top) the ordered 
c-phase is established. Note that the molecules do not 
form islands and the surface is thus completely covered 
in both cases. The overall photoemission intensity from 
the molecules is of course approximately a factor of two 
smaller than in the high coverage spectrum when normal- 
ized to the Ag background, simply because the number of 
molecules per area is smaller. Apart from that, the two 
spectra look identical: as in the case of NTCDA/Ag(lll), 
the former LUMO is partly occupied and appears be- 
low the Fermi level {Eb ~ 0.13 eV). Most important, 
for both samples the narrow peak at the Fermi level is 
clearly observed. The width (FWHM) of this signal is 
FWHM « 40 meV and is the same for both coverages. 
Again, the line width is strongly temperature dependent 
and increases following a linear behavior with A « 0.9 
when interpreted as due to electron-phonon coupling. 

Therefore it is evident, that long range lateral order 
of the adsorbate molecules is not necessary for the ex- 
istence of the narrow peak. Instead, it is the position 
of the LUMO which is important for the development of 
this additional feature. In the case of NTCDA/Ag(lll), 
there occurs a temperature and time dependent shift of 
the LUMO to higher binding energies, immediately con- 
nected to the disappearance of the narrow peak. 

With other words, the exact position of the LUMO is 
determined by hybridisation between the Bloch states of 
the substrate and the localized molecular states of the 
individual molecule. Such a scenario is described by the 
very versatile single impurity Anderson model (SIAM), 
which is based on the interplay between a localized mag- 
netic impurity, e.g. a transition metal or lanthanoid 
atom, with the conduction states of a metal [IS]. The 
SIAM has been applied successfully for the explanation 
of the physical properties of many transition metal and 
rare-earth compounds [50], and even for the quantita- 
tive estimation of many-body effects in adsorbate sys- 
tems (Newns- Anderson) [5T]. In particular, by use of 
numerical techniques [521 [S3| one can calculate the spec- 
tral function of such an impurity system, giving the so 
called Kondo resonance at the Fermi level. The max- 
imum position and the width of the Kondo resonance 
reflect the small energy scale fc^Tjf of the system, which 
determines all electronic and magnetic low-energy exci- 
tations in the system. Note that this model is based on 
a single impurity, coherence effects and the formation of 
an impurity lattice are not included and therefore not 
necessary for the appearance of the Kondo resonance. A 
typical temperature dependence for the SIAM line width 
derived from a non-crossing approximation (NCA) cal- 
culation [3S] is shown in fig. [2] {Tk='&^ K, off set by an 
additional broadening of 12 meV). A possible quantita- 
tive influence from coupling to phonons |16[ I17| is not 
discussed here. 

Indeed, the spectroscopic parallels between typi- 



cal Kondo (actually Heavy Fermion) systems, as e.g. 
CeCu2Si2 |38j or CeCue [21] are striking: a narrow fea- 
ture appears at the Fermi level with a strong temperature 
dependence in both line width and intensity. 

A quantitative determination of the Kondo scale ksTn 
from the temperature dependence requires a more de- 
tailed knowledge of the model parameters Coulomb cor- 
relation energy C7, single particle energy ej and hybridiza- 
tion strength V . However, a comparison with the spectra 
of inorganic systems allows a rough estimation of kgTi^. 
For example, the Kondo temperature must be signifi- 
cantly larger than for CeCue {Tk ~ 5 K), where the 
Kondo resonance vanishes nearly completely in the raw 
data for temperatures above 100 K. From the compar- 
ison with other Ce compounds, particularly CeSi2 j40j . 
we estimate Tk to be on the scale of roughly 100 K. This 
is in accordance with the low-T limit of the experimental 
line width (25 meV), which represents an upper limit of 
ksTx, and with the NCA simulation in fig. 2. 

The interpretation in the frame of the SIAM has fur- 
ther important implications, namely the meaning of the 
position of the LUMO. As shown in Fig. [T] the nar- 
row feature disappears, when the binding energy of the 
LUMO increases. Indeed, the investigation of other re- 
lated organic adsorbate systems [9] [HI [23l [29l [48] shows 
that only if the LUMO has its maximum below but close 
to the Fermi level — this might be seen as slightly above 
half filling — the narrow feature appears. This impor- 
tance is known from the 2D Hubbard model [SU [IS] with 
a band filling in the range of 0.6-0.7. 

Finally, for a "Kondo like" ground state, the correla- 
tion energy U must be large in comparison to the other 
model parameters, in particular to the hybridization be- 
tween conduction band and localized states V and to the 
single particle binding energy e, possibly given by the 
position of the LUMO maximum in this case, which is 
about 0.3 eV. An upper limit of 2 eV for U can be de- 
rived from the change in HOMO-LUMO separation from 
non-interacting molecules in multilayers to NTCDA at 
the interface. In addition, for the similar, yet larger, 
molecule PTCDA U is estimated to about 400 meV if 
adsorbed on Ag(lll) [56]. Following simple size argu- 
ments this value can be regarded as a lower limit for U 
of the smaller NTCDA. Moreover, a rough estimate of 
0.2 eV can be derived for V from the additional broad- 
ening that occurs for the LUMO in comparison to the 
HOMO (which shows weaker hybridisation [28] ) . 

In conclusion, we could demonstrate that photoemis- 
sion excitation spectra of chemisorbed organic molecules 
show a narrow feature at the Fermi level with a strong 
characteristic temperature dependence, that can not be 
explained by a mere electron-phonon coupling process. 
The appearance of the feature is strongly dependent on 
the position of the LUMO, i.e. on the band filling of 
the system, but does not depend on the long-range or- 
der of the monolayer. Therefore, we conclude that the 



observed spectral resonance at the Fermi level is a result 
of the coupling between the individual molecule and the 
conduction electrons of the metallic substrate. For an 
interpretation within the SIAM, the model parameters 
can be estimated to reasonable values. However, further 
quantitative studies by numerical methods are required 
to understand the many-body properties of these molec- 
ular systems. 
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